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Analogue models for CP violation in neutral-meson systems are studied in a general framework. 
No-go results are obtained for models in classical mechanics that are nondissipative or that involve 
one-dimensional oscillators. A complete emulation is shown to be possible for a two-dimensional os- 
cillator with rheonomic constraints, and an explicit example with spontaneous T and CPT violation 
is presented. The results have implications for analogue models with electrical circuits. 
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I. INTRODUCTION 



Many quantum systems display oscillatory behavior. 
Among the most interesting are the neutral mesons, 
where oscillations between particles and antiparticles can 
violate the product CP of charge conjugation (C) and 
parity (P) symmetries. In the K system, a small CP vi- 
olation is experimentally seen Q. It is associated with 
breaking of time-reversal symmetry T, with the product 
CPT being preserved Q). In fact, a complete formula- 
tion of CP-violating oscillations in the K system allows 
also CPT violation with T being preserved Hd]. A sim- 
ilar formulation can be developed for any of the neutral 
mesons K, D, B^, and B s ||. This more general situa- 
tion is of interest, for example, in the context of possible 
experimental signals from string theory || . 

Acquiring physical insight into the behavior of meson 
oscillations in the presence of T and CPT violation is 
worthwhile. One approach is to construct a simple ana- 
logue model in classical mechanics that displays the key 
features of meson oscillations. A priori, it seems most 
natural to adopt an intuitive picture based on an ana- 
logue model in which the meson and its antimeson are 
represented by two one-dimensional oscillators interact- 
ing through some weak coupling. Indeed, basic features 
of the CP-preserving physics can correctly be modeled 
in this way Q. However, modeling T violation is more 
subtle §,|. 

In this work, we investigate the issue of emulating both 
T and CPT violation in neutral-meson systems via mod- 
els in classical mechanics involving small oscillations. We 
obtain several no-go results, showing that complete em- 
ulation of the effective hamiltonian describing the time 
evolution of a neutral meson is impossible using mod- 
els with no damping or using models involving two one- 
dimensional oscillators with a large class of couplings. 
This confirms and extends earlier results of Rosner j^] . In 
contrast, models involving two-dimensional oscillations 
with appropriate constraints can display effects emulat- 
ing both T and CPT violation simultaneously. We give 
an explicit example in which the violation arises sponta- 
neously. 



In the next section, we present a few basic results 
needed for the subsequent analysis. Section III discusses 
some no-go results. The issue of spontaneous symmetry 
breaking is considered in section IV. The general analysis 
leading to a complete emulation of neutral-meson effec- 
tive hamiltonians, including an explicit model, is given 
in section V. Section VI summarizes the results and dis- 
cusses some open issues, including the implications of our 
results for the challenge of emulating CP violation with 
electric circuits. 



II. BASICS 

In this section, we introduce some basic results needed 
for the analysis in later sections. Following a discussion 
of relevant features of the neutral-meson systems, a few 
considerations appropriate for classical analogue models 
are presented. 

The four relevant neutral- meson systems are K°, D°, 
B%, and B®. In what follows, we denote by P° the strong- 
interaction eigenstate associated with any one of these. 
A general neutral-meson state is a linear combination 
of meson and antimeson wave functions. It can be repre- 
sented as a two-component object ^, with time evolution 
determined by a 2x2 effective hamiltonian A according 
to idt^f = A^ff. The eigenvectors Ps and Pl of A are 
the physical propagating states. The hamiltonian A is 
composed of a hermitian mass matrix M and a hermi- 
tian decay matrix T: A = M — \iT. Flavor oscillations 
and T violation are governed by the off-diagonal elements 
of A, while CPT violation is controlled by the difference 
between its diagonal elements. 

A widely used parametrization of A is |B[ 



A = 



-iD + E 3 Ex- iE 2 



Ei + iE 2 -iD - E 3 



(1) 



where D, E\, E 2 , E 3 are complex. In this parametriza- 
tion, T violation occurs when (EiE% — £7* £2) 7^ while 
CPT violation occurs when E3 ^ 0. In terms of real and 
imaginary components, T violation occurs when 



f 



(ReE 2 lmE 1 - ReE 1 ImE 2 ) ^ 0, (2) 
and CPT violation when either or both of 



Re E 3 ^ 0, Im E 3 ^ 0, 



(3) 



is satished. There are therefore three independent real 
quantities determining CP violation in neutral-meson 
systems. 

The basic goal in achieving the construction of a suit- 
able analogue model is to obtain an oscillating system 
in classical mechanics with a characteristic matrix repro- 
ducing the features of the effective hamiltonian A. It is 
therefore useful to consider the extent to which the form 
(|l|) of A can be modified without affecting the underlying 
meson physics. 

One flexibility in the form of A arises because the P° 
and P° wave functions are eigenstates of the strong inter- 
actions, which preserve strangeness, charm, and beauty. 
For a given system, the phases of the two wave functions 
can be rotated by equal and opposite amounts without 
observable consequences. This rotation induces a corre- 
sponding change in the phase of the off-diagonal compo- 
nents of A, which acts to mix E\ and E 2 but preserves 
the combination (E\E\ — E*E 2 ) measuring T violation. 
It follows that a satisfactory analogue model needs to 
contain features corresponding to this phase-independent 
measure of T violation. Note, however, that the phase 
rotation cannot mix components of the mass matrix M 
with components of the decay matrix T. 

Another flexibility relevant to the construction of an 
analogue model is the choice of basis for the meson sys- 
tem. The effective hamiltonian A in Eq. ([!]) is given in the 
P°-P° basis. However, other unitarily equivalent bases 
can also be chosen. A basis transformation by a unitary 
matrix U converts the state \f r into = and results 
in an effective hamiltonian A' = WAU. Appropriate 
choices for U can modify the location of the parameters 
for T and CPT violation in the effective hamiltonian, 
which may have some advantages in matching to an ana- 
logue model. Note that, like the phase-rotation freedom, 
the transformation by U cannot mix components of the 
mass matrix M with components of the decay matrix T. 

As an example, consider a CP-eigenstate basis Pi, P 2 
obtained from P , P° via the unitary transformation 



U = 



1 



V2 



1 1 



1 -1 



(4) 



In the new basis, the effective hamiltonian A becomes 
-iD + Ex E 3 + iE 2 



A = 



(5) 



E 3 - iE 2 -iD - Ei 



The elements of A are similar to those of A, except that 
Ei and E 3 have been interchanged. This conversion to A 



was used by Rosner and Slezak M to show that a mod- 
ified damped Foucault pendulum can be identified as an 
analogue model for T violation in the CP-eigenstate basis 
of a meson system. 

Other choices can be made. For example, combining 
the choice of CP-eigenstate basis with a phase rotation 
by exp(i7r/4) of the P° wave function and an opposite ro- 
tation of the P° wave function yields an effective hamil- 
tonian A given by 



A = 



-iD + E 2 E 3 - iEi 



E 3 + iEi -iD - E 2 



(6) 



This corresponds to a modification of A involving a cyclic 
permutation of the three parameters: (2?i, E 2 , E 3 ) — > 
(E 3 , Ei, E 2 ). 

In developing an analogue model, we adopt the notion 
that the behavior of the strong-interaction eigenstates 
P°, P° can be modeled classically by identifying them 
with harmonic oscillators in two generalized coordinates. 
The energies of the meson eigenstates are emulated by 
the oscillator frequencies, while the meson decay rates 
are paralleled by the oscillator dampings. 

Since the strong interactions preserve CPT, in the ab- 
sence of CP violation the two frequencies and decay rates 
are expected to be equal. The idea is to model the 
presence of CP violation by introducing appropriate cou- 
plings between the two classical oscillators. Following 
ideas concerning CPT violation in the context of con- 
ventional quantum field theory |5|,|l^], we regard it as 
desirable to obtain CPT violation spontaneously in an 
analogue model. In fact, we show below that it is also 
possible to generate T violation spontaneously. 

We limit attention to classical models involving small 
oscillations about equilibrium with linear equations of 
motion. Assuming harmonic behavior, the linear gen- 
eralized coordinates q\, q 2 can be combined in two- 
component form as Q — Re [Aexp(iwt)], where A is a 
complex two-component object. The equations of mo- 
tion can then be expressed by the action of a 2 x 2 matrix 
X(uj) on A, as XA = 0. The matrix X is the character- 
istic matrix of the classical oscillator. 

A suitable analogue model for a neutral-meson system 
is one for which the characteristic matrix X reproduces 
the features of the meson effective hamiltonian. In com- 
parisons between the analogue model and the meson sys- 
tem, it is useful to adopt a form for X analogous to that 
for A in Eq. (|l|) . We therefore introduce the parametriza- 
tion 



X = 



-iA + Br, B, - iB 9 



Bi+iB 2 —iA — B 3 



(7) 



where A, B\, B 2 , B 3 are complex. 
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The reader is cautioned that, despite the similarity of 
the parametrizations (Q) and ([?]), the detailed physical 
meanings of A and X differ. For instance, A involves a 
first-order time development while X involves a second- 
order one. A related point is that the meson state \& 
is intrinsically complex, with physical observables being 
related to the norm of the probability amplitude. In con- 
trast, the mechanical coordinate Q is real, and the cor- 
responding amplitude A is complex only as a convenient 
artifact. For example, opposite phase rotations between 
the two coordinates could produce a physically inequiva- 
lent result in the classical analogue model, whereas simi- 
lar phase rotations on the strong-interaction eigenstates 
have no physical effect in the meson system. 



III. MODELS WITHOUT DAMPING 

The intrinsic physical differences between the quantum 
system and the classical model might seem sufficiently se- 
vere to exclude emulation of subtle effects such as T and 
CPT violation. Indeed, several no-go results can be ob- 
tained concerning the existence of an acceptable analogue 
model for A under various circumstances. In this section, 
we discuss obstacles to the development of an analogue 
model in the absence of damping forces. The effects of 
dissipation are considered in section V. 

Consider first a lagrangian C describing small linear 
oscillations in a conservative classical-mechanical system. 
For present purposes, we restrict attention to a system 
with two degrees of freedom, although some of our for- 
malism and results apply more generally. 

Linearity implies that C is quadratic in the real gen- 
eralized coordinates Q(t) and the first time derivatives 
Q = dQ/dt. It can be written as 



X 



-u?T + iojG + V. 



(10) 



C = \Q T TQ + \Q T GQ - hQ T VQ, 



(8) 



where T, G, and V are square matrices of the same di- 
mension as Q. By inspection, T and V are symmetric, 
while G is antisymmetric. Since C is real, all three matri- 
ces can be taken real without loss of generality. We call T, 
G, and V the kinetic, gyroscopic, and potential matrices, 
respectively. Note that G violates classical time-reversal 
symmetry. 

The Euler-Lagrange equations of motion obtained from 
the lagrangian C are 



TQ + GQ + VQ = 0. 



(9) 



The gyroscopic matrix G does not represent damp- 
ing, despite its association with Q, because it is de- 
rived from a lagrangian and the corresponding general- 
ized force is conservative. For harmonic solutions with 
Q = Re [Aexp(iuit)] Eq. (||) becomes XA = 0, where the 
characteristic matrix X has the form 



This matrix is hermitian and so can be diagonalized 
with real eigenvalues. The normal- mode frequencies are 
obtained from the condition detX(w) = 0, which is a 
quadratic equation in to 2 . The absence of damping phys- 
ically implies that there are two real normal-mode fre- 
quencies, and this can be confirmed by inspection of the 
discriminant of the general solution for u> 2 . 

In terms of the parametrization Im) of X , we find: 



Re A 
ImA 



Im Bi = Im B 2 



ImB 3 = 



0, 



>2 

ReB 3 = -|w 2 (T la - T 22 ) + \{Vn - V 22 ), 
ReBi = -uj 2 T 12 + V 12 , 



ReB 2 



-LoG\ 2 



(11) 



The form of Eq. (11) permits several conclusions about 
the feasibility of constructing nondissipative analogue 
models for CP violation in neutral-meson systems. Next, 
we discuss these conclusions for T and CPT violation in 
turn. 

To begin, observe that Eq. ([llj) includes the result 
ImBi = lmB 2 = for all possible conservative classi- 
cal systems. In contrast, Eq. (||) implies that at least one 
of Im Ei and Im E 2 must be nonzero for T violation in a 
meson system. It follows by comparison of Eqs. (Q) and 
(0) that T violation in the P a -P° basis with the effective 
hamiltonian A cannot be emulated by any nondissipative 
classical model with two degrees of freedom. 

Transformation to some other basis for the meson wave 
functions offers more flexibility but remains insufficient. 
For example, in the CP-eigenstate basis a nonvanishing 
component Rei?i in the effective hamiltonian (|^) can 
be modeled with Re B3 , but no means to model Im E\ , 
ImE 2 , ImE^ exists. The point is that neither the phase- 
rotation flexibility nor the choice of wave-function ba- 
sis can mix contributions to the mass matrix M with 
those to the decay matrix T, as discussed in the previous 
section. Since Im-Ei, lmE 2 , ImE^ are contained in T 
while Re Ei, ReE 2 , Re S3 are contained in M, there is 
no means to convert one type of contribution to another. 

We conclude that it is impossible to emulate T viola- 
tion in a neutral-meson system with any nondissipative 
classical model having two degrees of freedom. In essence, 
a successful analogue model for T violation in a meson 
system must involve dissipation because T violation in 
the meson system itself intrinsically involves dissipative 
oscillations. 

The situation for CPT violation has both similarities 
and differences. Comparison of Eqs. ([!]) and (Q) shows 
that nonzero CPT violation in a meson system involving 
Re E3 can be emulated by a classical oscillator model for 
which ReB 3 ^ 0. The result (O) reveals that it suffices 
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to have a difference between the diagonal elements of ei- 
ther the kinetic or the potential matrix. This is straight- 
forward to achieve in a physical system. In contrast, an 
argument similar to that for T violation demonstrates 
that it is impossible to emulate CPT violation involving 
ImEs in a meson system with any nondissipative classi- 
cal model having two degrees of freedom. This can again 
be traced to the association of Im E 3 with the decay ma- 
trix r and hence with dissipation in the meson system. 

The strength of these no-go results suffices to show the 
need for dissipative classical oscillations. However, before 
turning to issues pertaining to spontaneous breaking and 
dissipation, we present some remarks about gyroscopic 
terms in the context of conservative systems. 

For a completely general emulation of neutral-mesons 
systems, we deem it desirable to construct an analogue 
model for which all eight real parameters in Eq. (0) are 
nonzero. The result (|ll|) shows that in the absence of 
damping a nonzero gyroscopic matrix G is needed to ob- 
tain a nontrivial Re£?2. In fact, this also holds in the 
presence of dissipation, as is shown in section V. Models 
without gyroscopic terms are therefore of lesser interest. 
However, gyroscopic terms appear in only a restricted 
class of models. In particular, there is no simple means 
of generating a nonzero G in models involving two cou- 
pled one-dimensional oscillators, as we discuss next. 

Prior to linearization, a general lagrangian describing 
two coupled one-dimensional oscillators involves a kinetic 
term for each oscillator and an interaction potential. No 
gyroscopic term is present. By definition, the kinetic en- 
ergy of a one-dimensional oscillator involves only one gen- 
eralized coordinate, so linearization of the kinetic pieces 
cannot generate the cross-coupling needed for a nonzero 
G. The potential term would therefore need to be the 
source of G. However, the gyroscopic term is linear in the 
generalized velocity, so any appropriate potential term 
must be velocity dependent. This leaves only a restricted 
class of possibilities. 

It can be shown that G makes no contribution to the 
hamiltonian, so any acceptable velocity-dependent po- 
tential must describe forces that do no work. Forces that 
do no work and are described by a velocity-dependent 
potential certainly exist. A standard example is the 
Lorentz force on a charged particle moving in a mag- 
netic field. One might, for example, consider a model in- 
volving two charged magnetic dipoles, each restricted to 
move along a one-dimensional curve so that the only pos- 
sible oscillations are indeed one-dimensional. The force 
i*2i oc vi x B2 on one dipole is determined by its velocity 
vi and by the field B2 of the other dipole, as needed. 
However, this fails to generate directly a nonzero G be- 
cause F21 ■ V\ = 0, so the force is orthogonal to the oscil- 
lation. 

In short, we find that it is difficult and perhaps impos- 
sible to emulate all eight parameters for a neutral-meson 
effective hamiltonian with any classical model involving 



two coupled one- dimensional oscillators. We conjecture 
that an impossibility proof could be constructed on the 
basis that G violates classical time-reversal invariance, 
which imposes severe constraints on one-dimensional sys- 
tems. In any event, it would be interesting to obtain an 
impossibility proof or to provide a simple counterexam- 
ple. 

The above result provides strong motivation to turn 
instead to an analogue model involving one two- 
dimensional oscillator. In this case, it is possible to gen- 
erate a nonzero G under suitable circumstances. 

Before linearization, the kinetic term of a two- 
dimensional model typically involves both generalized co- 
ordinates. If the equilibrium coordinates and configura- 
tion are independent of time (scleronomous constraints) 
and if there are no ignorable coordinates, then the ki- 
netic term is quadratic in generalized velocities |]llf and 
so no gyroscopic term emerges upon linearization. How- 
ever, for the special class of models with time-dependent 
(rhconomic) constraints, linearization of the kinetic term 
can generate a nonzero G matrix. For example, suppose 
the model involves small oscillations about a uniform mo- 
tion, characterized by a constant vq with dimensions of 
velocity. Then, linearization of the terms quadratic in 
generalized velocities can lead to expressions involving 
the product of v and the oscillation velocity Q. These 
are linear in Q and under suitable circumstances can yield 
a nonzero G. Indeed, the term 'gyroscopic' refers to the 
appearance of a nonzero G matrix in the description of 
small oscillations of uniformly rotating bodies. Note that 
the classical T violation necessary for gyroscopic terms 
emerges here as a result of the uniform motion. 

IV. SPONTANEOUS SYMMETRY BREAKING 

In this section, we discuss the issue of generating T and 
CPT violation spontaneously in the classical model. The 
idea is to seek an analogue model with an initial con- 
figuration displaying no T or CPT violation, but with 
a perturbative instability causing a natural dynamical 
evolution to a stable configuration in which small oscilla- 
tions violate both T and CPT. This parallels the mech- 
anism for spontaneous breaking of CPT in string field 
theory ||. Following the discussion in the previous sec- 
tion, we primarily restrict attention to the case of one 
two-dimensional oscillator without appreciable dissipa- 
tion. The situation for viscous damping is considered in 
the next section. 

Consider first a particle moving under the influence of 
gravity on the interior of a spherical bowl that rotates 
with constant angular speed Q about the vertical axis. 
The configuration with particle initially at the bottom of 
the bowl is a solution to the equations of motion. How- 
ever, an otherwise negligible friction between the parti- 
cle and the bowl makes this configuration perturbatively 
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unstable if SI 2 > g/a, where g is the gravitational accel- 
eration and a is the bowl radius. The position of stable 
equilibrium lies instead on the surface at a vertical dis- 
tance below the center of the bowl. This example, 
introduced by H. Lamb in his paper on kinetic stabil- 
ity in 1908 p2[ , provides a classical implementation of 
spontaneous breaking of rotational symmetry. 

By itself, this example is unsatisfactory as the basis for 
an analogue model for CP violation in neutral-meson sys- 
tems because no restoring force is associated with a small 
horizontal displacement from the equilibrium position on 
the bowl's surface. However, a more general surface with 
noncircular horizontal cross section can avoid this diffi- 
culty. One might, for example, consider a surface that 
is a spherical bowl at the bottom but that smoothly de- 
forms into a surface of uniform elliptical cross section as 
the height increases. In this case only two equilibrium 
points occur, located on the semi-major axis of the el- 
liptical cross section. In equilibrium, the particle rotates 
with the bowl. Small oscillations about either equilib- 
rium point are stable in both vertical and horizontal di- 
rections. 

For an explicit analysis in the case of a suitable general 
surface, we adopt cylindrical coordinates (p, <f>, z) with 
origin at the bottom of the bowl. Let the bowl's sur- 
face be determined by the equation /(p, (j>, z) — 0, where 
by assumption / satisfies all the necessary convexity and 
smoothness conditions. Then, the motion on the surface 
of a particle of mass m under gravity is determined by 
the lagrangian 



C = \m[p 2 + p 2 {^ + nf + z 2 ] 



mgz + mXf(p, (f>, z), 
(12) 



where A is a Lagrange multiplier. See Figure 1. 

For sufficiently large Q, spontaneous symmetry break- 
ing occurs. The equilibrium point (po, 4>q, Zq) is deter- 
mined by the equations / = 0, = 0, and poft 2 + 
gfp/fz = evaluated at the equilibrium point, where 
subscripts on / indicate partial derivatives. Taking p 
and po4> as generalized coordinates for small oscillations 
of frequency lu about the equilibrium point, a short cal- 
culation shows that the characteristic matrix has compo- 
nents 



X22 



-(1 + t 2 )uj 2 - n 2 

-X Q (f pp + 2Tf pz +T 2 f zz ), 
~2iVLuj - A o (/ P + r/^)/p , 



(13) 



where L = pofl 2 /g, Xq = g / f z , and the partial derivatives 
are again evaluated at the equilibrium point. Note the 
appearance of the off-diagonal gyroscopic terms ±2iflu, 
as expected. 




FIG. 1. Particle of mass m moving under gravity on the 
general surface / = of a bowl rotating at uniform speed fl. 

For suitable /, this characteristic matrix is sufficiently 
general to model all four parameters Imi, Re£?i, Re£>2, 
Re -B3. However, to generate a finite Re B\ in the absence 
of dissipation, either f pr f, or f$ z must be nonzero. 

A special case, used in the next section, is a bowl with 
horizontal cross sections near the equilibrium point form- 
ing ellipses of constant eccentricity e and semi-major axes 
with the same orientation. For definiteness, we consider 
the surface determined near the equilibrium point by 

f(p, 0, z) = p 2 (l - e 2 cos 2 0) - (1 - e 2 )(z/k) 2 ' n = 0. 

(14) 

This describes a bowl of uniform elliptical horizontal 
cross section and vertical cross section along the x axis 
determined by z = kx n . One of the two equilibrium 
points is at po = (fl 2 /nkg) 1 ^ 11 ' 2 ^, 4> = 0, z — kp^. 
Small oscillations about this point are described in the p 
and po4> coordinates by the characteristic matrix 



X 



-(1 + T 2 )lu 2 + (n-2)Q 2 



-2iQu 



-2iQu 



-oj 2 + e 2 n 2 



(15) 



The component X 2 2 involves a function e 2 ( 
for later purposes is defined generally as 6 



e 2 cos 2 



'0), which 
2 %) = 



4>q ). In the present case 

2\ 



e 2 (2cos 2 o - 1)/(1 

<f>o = 0, which gives e 2 = e 2 /(l — e 2 ). We note in passing 
that the oscillatory motions determined by X are stable 
for n > 2. 
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In this simple model, the term ReBi vanishes. The 
analysis in the next section shows that this can be 
avoided with the addition of appropriate dissipative 
terms. However, we note in passing that a nonzero Re B\ 
can be obtained without dissipation by a relatively sim- 
ple modification of the surface, involving a helical twist 
with height. The idea is to arrange matters so that the 
semi-major axis of the horizontal elliptical cross section 
rotates as z increases. It suffices to replace <f> in the bowl 
surface function / of Eq. ( |i"4| ) with a function tf> + 0(z). 
The equilibrium condition for 4> becomes </>o = —6{zo), 
and the characteristic matrix for small oscillations ac- 
quires an additional contribution. For example, choosing 
6{z) = tz/pq with constant r produces a characteristic 
matrix equal to the sum of Eq. (15) and a twist term X T , 
given by 



x T = f 2 rTn 2 



tT 1 



1 



(16) 



There is therefore a contribution to Re B\ determined by 
the twist constant r. 



V. MODELS WITH VISCOUS DAMPING 

In this section, we consider analogue models involving 
classical oscillators with dissipation. Since we are us- 
ing linear and homogeneous equations of motion and the 
correponding characteristic matrix to model the neutral- 
meson effective hamiltonians, we limit attention only to 
damping forces linear and homogeneous in the general- 
ized coordinates and velocities. We refer to such damping 
forces as viscous, although this is a somewhat broader 
definition than normally used by physicists. Note that 
dry friction can also give linear equations of motion, but 
typically leads to inhomogeneous terms and so is disre- 
garded here. 

The standard procedure in classical mechanics is to ob- 
tain viscous damping forces from a Rayleigh dissipation 
function, which is a symmetric quadratic form in the gen- 
eralized velocities. However, under special circumstances 
viscous damping can lead to linear homogeneous damp- 
ing forces involving also the generalized coordinates p3[ . 
This case is of direct interest in the present context. We 
therefore work here with a generalized dissipation func- 
tion T that can handle damping in a broader class of 
models 0. 

Up to irrelevant terms, we take J- to be a general 
quadratic expression in the small-oscillation variables Q 
and Q : 



T = \Q T RQ + Q T HQ. 



(17) 



As usual, the damping forces are determined by the 
derivative of T with respect to the generalized veloci- 
ties. The real symmetric matrix R contains the standard 



Rayleigh dissipation matrix for viscous damping, along 
with any contributions from other types of damping that 
generate forces linear in the generalized velocities. The 
real antisymmetric matrix H determines damping forces 
linear in the generalized coordinates. 

Combined with the Euler-Lagrange equations (^), the 
generalized dissipation function ( [l7| ) leads to equations of 
motion for small oscillations in the classical model given 
by 



TQ+{G + R)Q + (V + H)Q = 0. 



(18) 



A harmonic solution has the form Q = Re [Aexp(iujt)} as 
before, but in the presence of damping w is complex. We 
write u) = 2i:v + in = fx + in. In what follows, we also 
use to 2 = A 2 + 2i/iK, where A 2 = /i 2 — k 2 . To simplify 
the discussion, we take the magnitude of the damping to 
be sufficiently small that potential complications such as 
the issue of stability require no special attention. 
The characteristic matrix is 

X = -Tlj 2 + i{G + R)lu + V + H. (19) 

In terms of the parametrization in Eq. (Q), we find 

ReA — (jui(T u + T 22 ) - ± M (i? n + i? 22 ), 

W = -±A 2 (T n +T 22 ) 

-i K (i?ii+i? 22 ) + i(Un + U 22 ), 
ReBi = -A 2 Ti 2 - kR 12 + Via, 
ImBi = -2^kTi 2 + (J,Ri2, 
ReB 2 = -fJ,G%2, 
ImB 2 = —kG\ 2 + H\ 2 , 
RcB 3 = -iA 2 (T 11 -T 22 ) 

-\k{R 1x -R 22 ) + \{V 1x -V 22 ), 
JmB 3 = -iiK(T u - T 22 ) + ±M(E n - R 22 ). (20) 

Inspection of these expressions shows that a sufficiently 
general model can indeed emulate independently all eight 
real parameters in the effective hamiltonian for a neutral- 
meson system. Note that the parameter Re_B 2 is unaf- 
fected by dissipation, as mentioned in section III, imply- 
ing that a complete emulation of the neutral-meson sys- 
tem requires a nonzero gyroscopic term and therefore is 
most readily accomplished using a single two-dimensional 
oscillator. Note also that the damping force involving the 
matrix H contributes only to Im B 2 , whereas the matrix 
R affects all parameters other than B 2 . 

As an explicit realization of these ideas, we revisit the 
analogue model considered in section IV describing a par- 
ticle moving in a uniformly rotating bowl with surface 
function /. We suppose that the particle experiences an 
external viscous damping force. This might be imple- 
mented with a mesh bowl that allows resistance to the 
particle motion from the air or from some other static 
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fluid in which the bowl and particle are immersed. We 
take the generalized dissipation function for this resis- 
tance to be 



AX = -ATuj 2 + iARcu +AV + AH. 



(23) 



T = \mh[p 2 + p 2 (<p + n) 2 + i 2 ], 



(21) 



where the functional form of z = z(p, p, <f>, (f>) is under- 
stood to be determined from the bowl surface equation 
/ = 0. We also suppose that the parameter h, which con- 
trols the magnitude of the damping forces, is sufficiently 
small to avoid difficulties with stability. 

Inspection of the forces obtained from Eq. ( ^i"|) reveals 
that an additional constant damping force hp^fl in the <j> 
direction acts on the particle at equilibrium and moves 
the equilibrium position away from the previously de- 
termined location. For example, in the special case of 
a uniform elliptical horizontal cross section, the equilib- 
rium point is displaced from the apex of the ellipse. In 
general, the location of the new equilibrium point is de- 
termined by the simultaneous solution of the three equa- 
tions / = 0, f p + Tf z = 0, and / - Spo/z = 0, where 
£ = hT/Q. 

The dissipation function T in Eq. ( ^l|) describes the 
fluid resistance to the particle motion. For small oscil- 
lations, it includes both Rayleigh-type dissipation via a 
matrix R and damping linear in Q described by a matrix 
H. It thus implements the form of Eq. (|l7|). The as- 
sociated equations of motion can be derived, along with 
the accompanying characteristic matrix. We find that 
the components of the characteristic matrix are the sum 
of the corresponding components in Eq. ( [l3| ) with addi- 
tional terms given by the components of a matrix AX: 

AXu = ih(l + T 2 )uj, 

AX 12 = TScj 2 + A S(/ pz + r/„) - ihTZco, 

AX 21 = AX 12 + 2hn, 

AX 22 = -Y?u? + 2Ao£Wpo - X o^ 2 fz Z 

+ih(l + T l 2 )uj. (22) 



This result shows that the introduction of a relatively 
simple viscous damping force suffices to ensure that all 
four parameters Re A, Im£>i, I111B2, ImB^, can become 
nonzero. 

For the special case of the bowl with uniform el- 
liptical horizontal cross section described by Eq. (14), 
the incorporation of viscous damping via Eq. (^) re- 
sults in an equilibrium point at zq = f2 p^/ng, with 

Pi 

tan i 



= [(1 - e 2 )"Q 4 /(l - e 2 cos 2 o )"(nfcg) 2 ] 1/( "~ 2) and 
= —(1 — y/l — x)/a, where x = (1 — e 2 )a 2 and 
a = 2h/e 2 il. The requirement of real 4>o constrains the 
magnitude of h to \h\ < e 2 \Q\/2\ / l — e 2 . The corre- 
sponding characteristic matrix for small oscillations in 
the p and po4> coordinates is given by the sum of Eq. 
(15) with an additional matrix AX. In analogy with Eq. 
(19), AX can be taken to have the form 



Note that a putative term of the form AG is absent, as 
expected. The matrices AT, AR, AV, AH are: 



AT = 







-rs 



-it: 



Mi + r 2 ) 



-hTE 



AR = 



AV 



-hTj: /i(l + £ 2 



-(n- 1)MT 
(n-l)hn {n~2)h 2 



AH 




(24) 



(25) 



(26) 



(27) 



For small h and hence small £, the diagonal elements 
of the matrices AT, AV can be viewed as perturbations 
on the result (p"5|), which involves nonzero T, G, and 
V. However, the contributions from the off-diagonal ele- 
ments of AT, AV and from AR, AH are crucial for the 
complete emulation of a neutral-meson effective hamilto- 
nian. In particular, Eq. ( p0| ) shows that the parameters 
Re A, ImBi, lm.B 2 , lm.B% are all nonzero, as desired. 



VI. SUMMARY AND DISCUSSION 

This paper studied the emulation of indirect CP vio- 
lation in neutral-meson systems using oscillator models 
in classical mechanics. We obtained some no-go results 
for analogue models without damping and for ones in- 
volving two one-dimensional oscillators. The implemen- 
tation of spontaneous symmetry breaking was shown to 
be feasible. We proved that analogue models involv- 
ing one two-dimensional oscillator with rhconomic con- 
straints can suffice to emulate all eight real parameters 
in the meson effective hamiltonian, including the three 
describing physical T and CPT violation. 

We presented a specific analogue model that provides 
a complete emulation. It involves a particle moving un- 
der gravity on the surface of a uniformly rotating bowl 
of elliptical cross section in the presence of weak exter- 
nal viscous damping. The equations for small oscillations 
about an equilibrium point are determined by a charac- 
teristic matrix given as the sum of X in Eq. (|l5|) and AX 
in Eq. ( p3[ ) . The parametrization (Q) of this characteristic 
matrix, with parameters fixed by Eq. ( p0| ) , can be placed 
in one-to-one correspondence with the parametrization 
(El) in the P°-P° basis of the effective hamiltonian A for 
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a neutral- meson system. The correspondence is A <-> D, 
Bi «-> Bi, _B 2 «-> S2, B3 «-> £3. Correspondences also ex- 
ist with the effective hamiltonian (||) in the CP-eigenstate 
basis or, since the emulation is complete, in any other ba- 
sis. 

The results we have obtained leave open some interest- 
ing issues. One is the extent to which quantitative values 
of experimental observables in any neutral-meson system 
can be emulated in a realistic version of the models we 
have discussed. A satisfactory match would require re- 
producing the relative sizes of the values of the masses, 
lifetimes, and parameters for CP violation. Since the 
experimental data available on oscillations in the four 
neutral-meson systems range from being relatively com- 
plete for the K to limited for the B s , the degree of dif- 
ficulty in obtaining a satisfactory match varies consid- 
erably. In any case, the full flexibility of the analogue 
models is unnecessary because no CPT violation has been 
observed to date [p~5| P~7| - For the special case of the K 
system CP violation is observed to be small, so an em- 
ulation involving small dissipation is likely to be possi- 
ble. It would be interesting to determine the feasibility of 
constructing a quantitatively accurate model, including 
perhaps constructing a working prototype. 

A more ambitious task would be to explore the in- 
sights provided by the model about T and possible CPT 
violation with an eye to understanding its origin in na- 
ture and the neutral-meson systems. For example, it 
is intriguing that the no-go results strongly favor two- 
dimensional systems with rheonomic constraints. This 
suggests a preference for a dynamical origin of CP viola- 
tion. Similarly, it is interesting that the CPT violation 
in the analogue model emerges from a violation of rota- 
tion invariance. This would appear to correspond with 
the situation in conventional quantum field theory in the 
context of the standard model, where the known mecha- 
nism for CPT violation originates in the spontaneous vio- 
lation of Lorentz symmetry |Q and implies CPT signals 
in neutral-meson systems that depend on the orientation 
and magnitude of the meson momentum |[l9f . 

Another interesting topic, raised by Rosner M, is the 
emulation of T and CPT violation by electrical circuits. 
The general analysis we have provided in this work can 
offer some insights. A detailed analysis of this subject 
lies beyond the scope of this work, but in what follows 
we provide a few remarks. 

Suppose that each of the strong- interaction eigenstates 
P° and P° is modeled as an oscillating electric circuit, 
with CP violation regarded as a weak coupling between 
them. For definiteness, we view the two meson wave func- 
tions as corresponding to the charges qi(t), Q2 (t) flowing 
through the circuits as a function of time. As in the 
case of the analogue model in classical mechanics, the 
energies of the meson eigenstates are emulated by the 
oscillator frequencies while the meson decay rates cor- 
respond to the oscillator dampings. Inductances in the 



circuit replace masses in the mechanical model, inverse 
capacitances replace coupling constants in the potential, 
and resistances provide dissipation. 

The two-component meson wave function \1/ can be 
identified with a two-component object Q(t) formed from 
qi and q2, as for the case of a classical- mechanics model. 
We take the differential equations describing oscillations 
of the charges in the circuit to be linear in Q and its 
time derivatives Q — I and Q — I, where I is the two- 
component current. In the absence of dissipation, the 
equation governing the oscillatory behavior of Q is Eq. 
(0), where the matrices T, G, V are interpreted as charac- 
terizing appropriate properties of the circuit. This means 
that much of the analysis in section III applies in modi- 
fied form. In particular, the results obtained there reveal 
that a primary obstacle to a complete emulation of CP 
violation via electric circuits is the need for an antisym- 
metric matrix G coupling the currents in the two circuits. 

A suitable circuit realization of G requires a two-port 
device that is passive (no energy storage, increase, or dis- 
sipation). The antisymmetry implies that reciprocity is 
broken: a potential V applied across the first port would 
induce a current across the second differing in phase by 7r 
relative to the current induced across the first port when 
the same V is applied across the second. Remarkably, 
two-port devices of this type, called gyrators, have been 
the subject of some attention in the specialized electron- 
ics literature since their original invention by Tellegen in 
1948 pCfl . Moreover, a variety of network realizations of 
a gyrator exist pl| . 

We therefore suggest it is feasible to develop an elec- 
tric circuit emulating all eight parameters in the effective 
hamiltonian for a neutral-meson system, including both 
T and CPT violation. A gyrator would implement the 
crucial T-violating features of Eq. (^|) and in particu- 
lar a nonzero Rei?2 in Eq. (fill). As in the case of the 
classical-mechanics oscillators, it would also be necessary 
to include dissipation. This requires designing a circuit 
that incorporates suitable damping elements leading to 
the oscillatory behavior given by Eq. (|l8|). A general ma- 
trix R can be obtained by a suitable placement of resis- 
tors in the circuit, though producing a dissipation matrix 
of the form of H might be less straightforward. Devel- 
oping an electrical realization of spontaneous symmetry 
breaking would also be attractive. A circuit designed 
to exhibit all these features would make an impressive 
tabletop demonstration emulating T and CPT violation 
in neutral-meson systems. 

The results obtained in the present work may also have 
application in the emulation of other quantum oscilla- 
tions in physics. For example, it would be of interest 
to study analogue models for neutrino oscillations. A 
complete analysis for this case is likely to be more in- 
volved, partly because three neutrino species are known 
and the options for CP and CPT violation are corre- 
spondingly more complicated. Nonetheless, an explicit 
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analogue model in classical mechanics or with electric 
circuits could provide valuable insight. 
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